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A detailed theory of a phase diagram of a 2D surface superconductor in a parallel magnetic 
field is presented. A spin-orbital interaction of the Rashba type is knowrt^ to produce at a high 
magnetic field h (and in the absence of impurities) an inhomogeneous superconductive phase similar 
to the Larkin-Ovchinnikov-Fulde-Ferrel (LOFF) state with an order parameter A(r) oc cos(Qr). We 
consider the case of a strong Rashba interaction with the spin-orbital splitting much larger than 
the superconductive gap A, and show that at low temperatures T < 0.4T c o the LOFF-type state 
is separated from the usual homogeneous state by a first-order phase transition line. At higher 
temperatures another inhomogeneous state with A(r) oc exp(iQr) intervenes between the uniform 
BCS state and the LOFF-like state at gush ~ 1.5T c o. The modulation vector Q in both phases 
is of the order of g^iBh/vF- The superfluid density n v a v vanishes in the region around the second- 
order transition line between the BCS state and the new "helical" state. Non-magnetic impurities 
suppress both inhomogeneous states, and eliminate them completely at T c ot < O.lf. However, 
once an account is made of the next-order term over the small parameter cx/vf <s£ 1, a relatively 
long- wave helical modulation with Q ~ gfj,sha/v F is found to develop from the BCS state. This 
long-wave modulation is stable with respect to disorder. In addition, we predict that unusual vortex 
defects with a continuous core exist near the phase boundary between the helical and the LOFF-like 
states. In particular, in the LOFF-like state these defects may carry a half-integer flux. 

PACS numbers: 74.20.Rp, 74.25. Dw 



I. INTRODUCTION 

There are experimental indications^ in favor of exis- 
tence of superconductive states localized on a surface of 
non-superconductive (or even insulating) bulk material. 
Such a superconductive state should possess a number 
of unusual features due to the absence of the symme- 
try "up" v/s "down" near the surface: the condensate 
wave- function is neither singlet nor triplet, but a mix- 
ture of both^; the Pauli susceptibility is enhanced at 
low temperatures^ (as compared with the usual super- 
conductors); the paramagnetic breakdown of the super- 
conductivity in a parallel magnetic field is moved to- 
wards much higher field values due to a formation of 
an inhomogeneous superconductive stated similar to the 
one predicted by Larkin-Ovchinnikov and Fulde-Ferrel^ 
(LOFF) for a ferromagnetic superconductor. All these 
features steam from the chiral subband splitting of the 
free electron spectrum at the surface, due to the pres- 
ence of the spin-orbital Rashba term 7 -; the magnitude of 
this splitting apF is small compared to the Fermi energy 
but can be rather large with respect to other energies in 
the problem. The line of transition from normal to (any 
of) superconductive state T c (h) was determined in j4j; 
however, the nature of the transition between the usual 
homogeneous (BCS) superconductive state at low fields 
and the LOFF-like state at high fields was not studied. 

In this paper we provide a detailed study of the phase 
diagram of a surface superconductor in a parallel mag- 
netic field h (a brief account of our results was pub- 



lished in Ref. Q). We show that at moderate values of 
h ~ T c j \ib the behavior of this system is rather different 
from 2D LOFF model of @. Namely, we demonstrate 
the existence of a short-wavelength helical state with an 
order parameter A cx exp(iQr) (where Q _L h and 
Q ~ fJ-Bh/vp) in a considerable part of the phase dia- 
gram, which is summarized in Fig. [TJ The line CT is the 
second-order transition line separating the helical state 
from the homogeneous superconductor. Below the T 
point a direct first-order transition between the homoge- 
neous and the LOFF-like state takes place. The line TO 
shown in Fig. [1] marks the border of metastability of the 
BCS state; we expect that the actual first-order transi- 
tion line is (slightly) shifted towards lower H values. The 
line ST' marks the second-order transition between the 
helical and the LOFF-like state. The above results are 
valid within the leading (the zeroth-order) approximation 
over the small parameter a/vF <C 1; an account of the 
terms linear in a/vp <C 1 leads, in agreement with (To| . to 
the transformation of the uniform BCS state into a long- 
wavelength helical state with a wave vector q ~ aH/v F 
at the lowest magnetic fields. Therefore, the CT line is 
actually a line of a sharp crossover (with a relative width 
of the order of cy/vf) between the long- and the short- 
wavelength phases. 

The rest of the paper is organized as follows. In Sec.lITl 
we introduce a model of a spin-orbital superconductor in 
a parallel magnetic field, with a hierarchy of the energies 
£p » apF S> T c . In Sec. Mil we derive the Ginzburg- 
Landau functional for an inhomogeneous ground state. 
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On the T C (H) line we find two critical points: the Lif- 
shitz point £ and the symmetric point 5, and demon- 
strate the existence of a "helical" state with an order 
parameter A(r) oc exp(iQr) and a large Q ~ H/vp. The 
point S on the T C (H) line is special in the sense that 
there the order parameter symmetry is enhanced to U(2) 
from the usual U{\). That leads to unusual vortex tex- 
tures discussed in Sec. IIVI In particular, vortices with 
half-integer flux are predicted for the LOFF-like state. 
In Sec. IVIwe derive the two stationary conditions for the 
helical state, which determine the equilibrium A and Q. 
The latter allows to establish the boundaries of stability 
of the BCS and helical state: the Lifshitz line terminat- 
ing in the critical Landau point T, and a line starting in 
the symmetric point S. We show that the helical state 
and the parity-even (stripe) phase are separated by two 
phase transitions of the second order and an intervening 
additional superconducting phase. In Sec. IVII we prove 
that the electric current is zero in the ground state; then 
we show that the supercurrent response to the vector po- 
tential component A y = AQ/Q vanishes at the CT line, 
whereas within the helical state both components of the 
superfluid density tensor n s are of the same magnitude 
as the superfluid density in the BCS state. In Sec. IVIII 
we explore the influence upon the phase diagram of the 
terms linear in q/dj? -C 1. We show that at low mag- 
netic fields the ground state is realized as a weakly helical 
state with zero current. A special geometry is proposed 
when an oscillating supercurrent in the ground state of 
the helical phase may flow. In Sec. IVIIII we study the 
effect of the non-magnetic impurities on the phase dia- 
gram. We show that in the relatively clean case the para- 
magnetic critical field is quickly suppressed by disorder 
and the position of the Lifshitz point is shifted towards 
higher magnetic field values. We find the critical strength 
of disorder above which all short-wavelength inhomoge- 
neous states are eliminated from the phase diagram; in 
terms of an elastic scattering time r this condition reads 
t c = O.llh/T c o. At t < t c the only phase which sur- 
vives is the "weakly helical" state with Q — 4aH/v F ; in 
this regime the paramagnetic critical field starts to in- 
crease with disorder, and at r -C t c we find H c oc r -1 / 2 . 
In Sec. IIXI we go beyond the mean-field approximation 
and study the modifications of the transition line T C (H) 
due to the Berezinsky-Kosterlitz-Thouless vortex depair- 
ing transition. We demonstrate that vortex fluctuations 
are strongly enhanced near the points C and S, leading 
to local downward deformations of the actual Tbkt(H) 
line. 



II. MODEL OF A SPIN-ORBITAL 
SUPERCONDUCTOR 

Near the surface of a crystal translational symmetry 
is reduced and inversion symmetry is broken even if it is 
present in the bulk. (The component of the electron mo- 
mentum p parallel to the interface is conserved because 



of the remaining 2D translational symmetry.) As a re- 
sult a transverse electrical field appears near the surface. 
The electron spin couples to this electric field due to the 
Rashba spin-orbit interaction 7 



a a x p 



(1) 



where a > is the spin-orbit coupling constant, n is a 
unit vector perpendicular to the surface, a = (a x , <J y ,cr z ) 
are the Pauli matrices. This interaction explicitly vi- 
olates inversion symmetry. The electron spin operator 
does not commute with the Rashba term, thus the spin 
projection is not a good quantum number. On the other 
hand, the chirality operator [a x e\ • ft commutes with 
the Hamiltonian. Here e — p/p is the momentum di- 
rection operator with eigenvalues e p = (cos tp p , sinyj p ), 
where <p p is the angle between the momentum of the 
electron and the £-axis. The chirality operator eigen- 
values A = ±1 together with the momentum constitute 
the quantum numbers of the electron state (p, A). The 
Rashba term (JTJ) preserves the Kramers degeneracy of the 
electron states, thus the states (p, A) and (— p, A) belong 
to the same energy. 

In this paper we consider the simplest model 3 of a sur- 
face superconductor: a BCS model for a two-dimensional 
metal with the Rashba term (JTJ) , in the limit apF ^> T c . 
The Hamiltonian written in the coordinate representa- 
tion reads 

H = J yj+(r)hM^d 2 r^^ J ^+^ a (ff, (2) 
with the one-particle Hamiltonian operator 
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X P 



ft- gfiBh ■ c? a p/2 , (3) 



where m is the electron mass, a, (3 are the spin indices 
and P — — iV — -A(r) is the momentum operator in 
the presence of an infinitesimal in-plane vector poten- 
tial A = A(r), e < is the electron charge. We have 
included into the Hamiltonian the Zeeman interaction 
with a uniform external magnetic field h parallel to the 
interface, assuming h to be in the x-direction. The vec- 
tor potential of such a field can be chosen to have only 
the z-component, therefore it decouples from the 2D ki- 
netic energy term. fiB is the Bohr magneton and g is the 
Lande factor. Hereafter we use a notation H = g^iBh/2. 
The electron operator can be expanded in the basis 



of plane waves: ip a (r) = J2 P \ eipTa ap, and the one- 
particle part of the Hamiltonian ([2]) in the momentum 
representation can be written as a sum of Hq and H e m ■ 



H r , 



2m 



-1 



fa/3 



3 A a^. 



H ■ (Tap I a/3 p 



(4) 



3 



17 T 




in-plane magnetic field 



FIG. 1: A phase diagram that shows: a superconducting phase transition line T C (H) (solid) and two second-order phase 
transition lines in the clean case, an CT line between the homogeneous (BCS) and the helical (h.) state, and an ST' line of 
stability of the helical state. The dotted line going downwards from point T to point O marks the absolute limit of stability of 
the BCS state. The cross indicates the Lifshitz point £ and the circle indicates the symmetric point S. The line of transition 
into the gapless superconductivity H — A is marked with a dashed line. 



Here the current operator is 

j = —e(p/m — a[a x ft]) — 



2mc 



(5) 



The Hamiltonian Hq can be diagonalized by the trans- 
V\a(p)o,\p with the two-component 

(6) 



formation a D 
spinor 



^ (P) " V2 (iAexp(iV P (ff)) 



where 



ip p (H) = arcsin 



ap y — H 



^/{ap) 2 - 2ap y H + H 2 



(7) 



The eigenvalues of the Hamiltonian (j4|) , corresponding to 
the chiralities A = ±1, are 



ca(p) - P 2 /2m - {ap) 2 - 2a Py H + H 2 , (8) 

thus at H = the equal-momentum electron states are 
split in energy by 2apF- The two Fermi circles, corre- 
sponding to the different chiralities, have Fermi-momenta 
Pp = y/2myL + m 2 a 2 ima, where \i 3> ma 2 is the chem- 
ical potential. The density of states on the two Fermi cir- 
cles is almost the same, v+ = I 1 ± — ] . In the main 




FIG. 2: When an external magnetic field H <C apF is applied 
in the ^-direction, the two Fermi circles corresponding to the 
different chiralities A = ±1 are shifted in opposite y-directions 
by a small momentum Q = ±H/vf- The circle corresponding 
to one of the chiralities is larger than the other by the relative 
amount of o/vf- 



part of the paper we neglect the difference v + — V- . The 
effects related to v + ^ v_ will be considered in Sec. IVIII 
When an external magnetic field is applied, these two 
Fermi circles are displaced in opposite y-directions by a 
momentum Q = ±H/vp, as shown in Fig. [5] 

The two-particle pairing interaction in Hamiltonian ([2]) 
in the momentum representation reads 



-E 

pp'q 



a ap+q/2 a /3-p+q/2 a /3-p'+q/2 a ap'+q/2' ( 9 ) 
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and can be simplified in the chiral basis ((6|), assuming 
H <C apF ^ /•*■ I n the long- wavelength limit q <C Pf h 
can be factorized as 



ff int = -^i+(q)i(q), 



(10) 



where the pair annihilation operator 



i(q) - £ Ae4VPa A-p+q/2 a Ap+q/2- (11) 

p,A 

Here Ae I¥3p is the wave function of the Cooper pair in the 
chiral basis and it changes sign under the substitution 
— p for p. 

To calculate the thermodynamic potential Q = 
—ThiZ, we employ the imaginary-time functional inte- 
gration technique with the Grassmanian electron fields 
OAp, «Ap and introduce an auxiliary complex field A(r, r) 
to decouple the pairing term Hi nt , cf. ll]. The resulting 
effective Lagrangian is 

L[a, a, A, A*] = ^ «a p (~d T - e A (p)) a Xp + 



p,A 



q p,a 



+ A* Ae 4Vp OA -p+q/2aA : p+q/2 



a A,p+q/2aA,-p+q/2 

(12) 



The thermodynamic potential = —TlnZ describes a 
system in equilibrium, where Z is the grand partition 
function. We express O as a zero-field limit of a generat- 
ing functional £1 = Cl[r], fj]\ 



exp 



°^^=/^A*expf-^^ A ' A * 



T 



T 



-JvaVaD AVA* exp (J ' [L[a, a, A, A*] - 
+ /^2(v(Tp)a(Tp) + o(rp)?j(Tp))] dr 



(13) 

Below we will work within the mean-field approxima- 
tion, which is controlled by the smallness of the Ginzburg 
number Gi. However, for a clean 2D superconductor 
Gi ~ T c /Ep may be non-negligible (we discuss the fluc- 
tuational effects in the end of this paper) . The mean-field 
approximation is equivalent to the saddle-point approx- 
imation for the functional integral over A and A* de- 
fined in the first line in Eq. (|13[) . In other terms, we 
will study the minima of the functional fi[A, A*], which 
comes about after integration over the Grassmanian fields 
in the functional integral defined in the second line of 

Eq. m 



5Cl[A, A* 
SA(rr) 
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FIG. 3: A Cooper loop with transferred momentum Q. 

To evaluate the thermodynamic potential 0[A, A*], we 
will use the Green's function method. The electron 
Green's function is defined as a variational derivative of 
the generating functional: 



G{tt,t't') = 



5r)(rr)dri(T'r') 



(15) 



In the next Section we determine the line of the su- 
perconducting transition T C (H), and locate two special 
points on this line, C and T, which designate the bound- 
aries of different superconductive states. 



III. SUPERCONDUCTING PHASE 
TRANSITION 

Near the phase transition from a normal metal to a su- 
perconductor the order parameter A(r) is small. There- 
fore the thermodynamic potential f2 may be expanded 
in powers of A(r) and its gradients. This is known as 
the Ginzburg-Landau functional. It has been shown by 
Barzykin and Gor'kov 4 , that the ground state can be in- 
homogeneous in the direction perpendicular to the mag- 
netic field. We consider the order parameter as a super- 
position of a finite number of harmonics: 



A(r) =^A Q( (r)exp (iQ t f) . 



(16) 



where Aq { (r) are slowly varying envelope functions. The 
corresponding Ginzburg-Landau functional is 



,Q 3 Q k Qi Aq; (r)^*Q J (r)A Qk (r)A* Ql (r) 

ijkl 

xS Qi+Qk-Q s -Qi + 



i 



-in^--^A x {r) j A Q ,(f) 



d 2 7 



(17) 



The coefficient oqq is given by the Cooper loop di- 
agram with transferred momentum Q. The coefficients 
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Pqqqq an d Pqq-q-q are determined by four-Green's- 
f unction loop integrals: 

&QQ = Jj~^ E G a P + Q/2) G\ (-w, -p + Q/2) , 



cj,p,A 



(18) 



/3QQQQ = T J2 G l (w. P + Q/2) G a (-w, -P + Q/2) , 

0QQ-Q-Q =TJ2GI(UJ,P + Q/2) X 

LJ,p,A 

xGa (-w, -P + Q/2) G\ (-w, -p + 3Q/2) , 



1 92 



-00 - 2 cV 



(19) 



(20) 



Here the normal state Green's function in an external 
in-plane magnetic field H <C ap F is 

1 



G\ (w,p) 



iu> — £\ (p) — \H sin y> p ' 



(21) 



where u> — 2n{n + |)T is the Matsubara frequency, n is 
an integer, and the quasiparticle dispersion 



£a(p) = p 2 /2m - XapF - H 



(22) 



is assumed to be small compared to cv^f- The integrals 
over the momenta in (fT5|) and (fTO)) are calculated in the 
semiclassical approximation: 



d 2 p 

(27? 



^a(ef) 



00 />27T 1 



In this Section we will calculate all diagrams in the zeroth 
order over the parameter a/v F <C 1 (i. e. we approxi- 
mate v\{e F ) ~ K £ f)); the effect of the terms of the order 
of 0(a/v F ) will be discussed in Sec. IVIII below. 

Integrating over the momenta p in the Ginzburg- 
Landau functional coefficients gives 



aQQ 



- - 7rv(e F )T 



(24) 



Pqqqq — 



4 

u(e F ) 



(w 2 + fff) 5 / 2 ' 



w>0,A 



(2^ 2 + ff 2 ) tf A 



2 '" ^ u 2 (u; 2 + H 2 ) 3 / 2 v F Q' 

w>0,A A 



(25) 
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FIG. 4: Diagrams corresponding to the terms of the fourth 
order in A in the Ginzburg-Landau expansion: /3qqqq|Aq| 4 
and (3qq-q-q\Aq\ 2 \A~q\ 2 . 
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FIG. 5: Cooper pairing wave vector Q. The circle indicates 
the symmetric point 5*. 



where 



H X =XH + v F Q/2. 



(27) 



Note, that in the Ginzburg-Landau functional (fTT]) the 
coefficients (3qqqq = (3-q-q-q-q correspond to the 
terms |Aq| 4 and |A_q| 4 ; and the coefficient correspond- 
ing to the term |Aq| 2 |A_q| 2 is a sum of the four equal 
coefficients Pqq-q-q = 0q-q-qq = P-q-qqq = 
P-QQQ-Q- 

The condition cvqq = determines the second-order 
transition line (if (3qqqq > 0) between the normal metal 
and the superconductor: 



— = v(e F )T m&x 
U Q 



uj>0,\ 



^uj 2 + (\H + v F Q/2) 2 



(28) 



Equation [[28]) determines the shape of the phase tran- 
sition line T C (H) between the normal and superconduct- 
ing states. Depending upon H, the maximum over Q 
in the r.h.s. of Eq. (|28[) is attained either at Q = 
or on nonzero Q = ±\Q\- The position of the T C (H) 
line found via numerical solution of Eq. I|28p is shown 
in Fig. [TJ where both temperature T c and in-plane mag- 
netic field H are normalized by the critical temperature 
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at zero magnetic field: T c q = 2u>d exp(— Y/vU + 7)/7r, 
where 7 = 0.577 is the Euler constant. The line T C (H) 
recovers two asymptotics found in 0]: 



log 



T C (H) 7a3)H 2 



T C {H) _ 
T, 



^ 2 T 2 



in the limit H/T t 



cO 



0. 



ttT, 



cO 



in the limit H/T c0 -> 00. (29) 



Whereas at low H the superconductive solution is uni- 
form, Q = 0, in the high field limit one finds^ Q = 
2H/vf- The Lifshitz point C separates Q = and Q ^ 
solutions on the T C (H) line and is the end of the second- 
order phase transition line between the two supercon- 
ducting phases. In order to determine the position of the 
C point, we note that cvqq, cf. Eq. (124[) . is symmetric 
under the change — Q — > Q, and thus it has always an ex- 
tremum at Q — 0. Therefore the position of the Lifshitz 
point C should satisfy the equation 



d 2 a QQ 



dQ 2 



= 0. 



Q=0 



Numerical solution of the above equation gives the value 
(H L ,T L ) = (1.536, 0.651)T c0 , with the ratio H L /T L w 
2.36. 

Fig. [5] shows the Cooper pairing wave vector Q on the 
T C (H) line as a function of the in-plane magnetic field H . 
Near the C point the wave-vector Q contains a square- 
root singularity vfQ{H) ~ y/ H 2 — H 2 , typical for the 
behavior of the order parameter near a second-order tran- 
sition. In the high-field limit H/T C (Q) — > 00 the behavior 
of the wave-vector Q is given by the asymptotic expres- 
sion 



v F Q = 2H 



7C(3)e 2 ''iJ 3 ' 



(30) 



Note, that Q — 2H/vp is the momentum shift of the two 
A = ±1 Fermi circles in a parallel magnetic field. 

Near the T C {H) line the coefficient a can be approxi- 
mated as 



a(T, H) w v(e F ) 



where Y(T, A) 



T - T C {H) 



2T 



Y,Y{T C {H),H X ), (31) 



1 



TJ2 



■^p J d£ sech 2 -^ 2T A ^ s ^ ne Yosida function, 
uj n = irT(2n + 1) is the Matsubara frequency. At 
H > Hl an inhomogeneous superconductive phase is 
formed below the T C (H) line. Eq. I[28p determines the 
absolute value of the equilibrium wave vector \Q\ (the 
direction of Q is perpendicular to H), therefore two 
harmonics may contribute to A(r) just below the T C (H) 
line: A(y) = A+e lQy + A_ e - iC ^. Below T C (H), the 
density of the thermodynamic potential f2 is lower in 
the superconductive state than in the normal one by the 
amount 



Q sn = a(T,H)\A\ 2 + 
&(T,tf)|A| 4 +/? (T,if)(|A H 




in-plane magnetic field 



FIG. 6: The ratio of (3qq-q-q and Pqqqq coefficients. The 
symmetric point <S, marked on the figure with a circle, is 
defined as a point where the ratio attains the value 1/2. The 
Lifshitz point C is the point where the ratio equals 1, since 
Q = there. 



where |A| 2 = |A + | 2 + |A_| 2 . Eq. ([32]) was obtained from 
Eq. (TP?)) by keeping of only two harmonics: |A + | 2 = 
AqAJj, |A_| 2 = A_ Q A1 Q . Comparing Eqs. QU) and 
([521 gives 



Ps(T,H) = -Pqqqq + Pqq-q-q, 



(33) 



Pa(T,H) — -Pqqqq - Pqq-q-q, 



where Pqqqq and Pqq-q-q are the four-Green's- 
function loop integrals (|25p . In the symmetric point iS, 
where p a {T c {H),H) = 0, the free energy ([3"2"| depends 
upon |A| 2 only, and thus is invariant under U(2) rota- 
tions of the order parameter spinor (A + , A_). The coor- 
dinates of the S point are: (H S ,T S ) = (1.779, 0.525)T cO , 
the corresponding wave vector is vfQs = 2.647T C (0). 
At H < Hs we find p a < 0, and the free energy at 
T < T C (H) is minimized by the choice of cither A + = 
or A_ = 0, both corresponding to the helical state. 
At H > Hs Pa > and the free energy minimum is 
achieved at |A+| = |A_|, i. e. the LOFF-likc phase with 
A(y) oc cos(Qy) is the stable one at high field values. 



IV. UNUSUAL VORTEX SOLUTIONS NEAR 
THE SYMMETRIC POINT S. 

A. General considerations: extended U(2) 
symmetry and vortices 

In this Section we discuss the peculiar properties of 
the superconductive vortices, which appear due to the 
extended symmetry of the order parameter near the sym- 
metric point S of the phase diagram. The free energy of 
A_ | 2 ) 2 , (32) the superconductor in the vicinity of the symmetric point 
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is given by the Ginzburg-Landau functional 



d 2 r[ 



A: A, 



(ci (idi + —i 
/U|A+| 

+/3 s (Ag-|A + | 2 -|A_| 2 ) 2 -^),(34) 



,2 | A _|2 )2 _ 



idi + yA: ! A _ 



where i — x,y, and A 2 , = —a/2(3 s is the equilibrium 
value of the order parameter. The fourth-order term in 
the Ginzburg-Landau expansion (|32|) can be divided into 



a symmetric and an anisotropy part /3 a (T, H)(\A + \ 2 — 
|A_| 2 ) 2 , and in the symmetric point the coefficient 
Pa{Ts, Hs) = 0. The coefficients c x , c y are given by the 
expressions (|26[) . and their ratio in the S point is equal 
to c x /c y = 1.72. Eq. (|34[) can be written in a rotation- 
ally symmetric form via an area-conserving transforma- 
tion by stretching and contracting the two coordinates 
x — > x(c x /c y ) 1/4: and y — > y(c y /c x ) 1/A . 

The minimum of the free energy (|34[) in the symmet- 
ric point (where /3 a = 0) is achieved in the homoge- 
neous state under the condition |A + | 2 + |A_| 2 = A 2 . 
In normalized variables Z\ = A + /A and 22 = A_/A 
the order parameter spinor (z±, z-i) spans the sphere S 3 : 
\zi 1 2 + |«2 1 2 = 1; an d is equivalent to a four-component 
unit vector TV. This normalization allows us to write 
the gradient part of the free energy 1|34[) as a non-linear 
sigma-model 



grad 



(35) 



where p s — ^-y/CxCy is defined through coefficients of 
the Ginzburg-Landau functional (|3~4"|) . 

Precisely at the symmetric point, the gradient func- 
tional ([33)) governs the system's behavior at length scales 
L larger than the order parameter correlation length 
£(T). At (3 a 7^ its applicability is limited from the 
large scales also. Namely, L should be smaller than the 
temperature-dependent "anisotropy length" 



L an (T)=C(T) 1 



A 

I A, 



(36) 



which is determined by comparison of the gradient term 
and the anisotropy term in the full free energy func- 
tional (34}. 

On the left of the symmetric point (/3 a < 0) the min- 
imum of the energy is achieved at either \z\\ = 1 or 
\z%\ — 1, leading to the degeneracy manifold of the order 
parameter, S 1 ®Zi- On the right of the symmetric point 
(f3 a > 0) the degeneracy manifold of the order parameter 
is S 1 (£> S 1 . In order for the gradient energy of a physical 
defect to be finite, at large distances from the defect core 
the order parameter should belong to the corresponding 
degeneracy manifold. On the other hand, at relatively 



small distances r < L ani the whole S 3 sphere is available 
for the order parameter configurations. 

The two-dimensional x-space is topologically equiva- 
lent to a sphere S 2 with a boundary at infinity S . A 
physical defect is described as a mapping of a disk on the 
real plane R 2 with a boundary S 1 (which encloses the 
defect) on the degeneracy manifold of the order parame- 
ter S 3 . Due to boundary conditions at infinity (imposed 
by the anisotropy) the mapping S 2 — > S 3 is accompa- 
nied either by the mapping S 1 — > S 1 (at (3 a < 0, i. 
c. on the left from the S point), or by S 1 — > S 1 <g> S 1 
(at f3 a > 0). Therefore the topological defects are de- 
termined by the nontrivial elements of the homotopy 
group ^(S* 3 ,^ 1 ) = Z (on the left from the S point) 
or n 2 {S 3 ,S 1 <£> S 1 ) = Z <£> Z (on the right from the S 
point). Note that in the absence of any anisotropy there 
would be no stable topological defects since TT2(S 3 ) = 0, 
i. e. any configuration of the order parameter could 
be transformed into a homogeneous state. The general 
approach to the classification of vortices with a nonsin- 
gular core by means of the relative homotopy groups 
TT2(R,R), described above, was first introduced by Mi- 
neev and Volovik [T3|; a review of the approach can be 
found in (T3 |. Some explicit solutions for nonsingular 
vortices are presented in the review [lil ]. 

For an actual calculation it is convenient to employ the 
Hopf projection, which splits the order parameter spinor 
z G S 3 , parameterized here as 



-iip/2, 
up/2 



is 0/2 \ 
,6/2 ) 



(37) 



into an N = 3 - unit vector n = z^uz € S 2 sphere, 



sin 6 cos ip 
sin 9 sin ip 
cos 9 



(38) 



parameterized through the Euler angles on the S 2 sphere, 
and a total phase \ € 1/(1), canonically conjugated to the 
electron charge. 

Each configuration n(x) defines a mapping of the co- 
ordinate plane (equivalent to S 2 ) on a sphere n 2 = 1, i. 
e. the mapping S 2 — > S 2 . These mappings are charac- 
terized by an integer "topological charge" 



Q = -j— { n \d x n x d y n] d 2 r, 

47T J 



(39) 



which is related with the circulation of the vector 



= -iiz^d^z - zd^J) = (d^x - -d^cosO). (40) 
Indeed the following identity can be proven: 

Q = lfA.dl-lfA.dl^§--C , (41) 

Coo Co 
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where Coo is the closed loop at infinity, and Co is the in- 
finitesimal closed loop just around the vortex singularity 
point. The last equality in Eq. (j4Tj) relates the topo- 
logical charge Q and the magnetic flux connected with 
the vortex defect (in units of the superconductive flux 
quantum <E>o). 

Now the gradient energy can be represented as a 
sum of the gradient energy of the n-field and the kinetic 
term: 



large energy 



d 2 f. 



(42) 



Below we consider separately vortex solutions in the he- 
lical state realized at (3 a < and in the LOFF-like state 
at p a > 0. 



B. Non-singular vortices in the helical state. 

In the region on the left of the symmetric point {(3 a < 
0) the energy minimum in the bulk of the film is attained 
at either Z\ = 1, Z2 = or vice versa; we choose the first 
solution for further discussions. Then the order parame- 
ter is proportional to e^ - ^! 2 *) and an elementary vortex 
corresponds to the rotation by 2ir of the "effective" pa- 
rameter's phase x~f/^ = <Pi where <p is the azimuthal an- 
gle on the plane. Near the vortex core, however, one can 
construct solutions which contain both Z\ and Z2 com- 
ponents. Now we show that such a solution has a lower 
energy than that of a standard singular vortex (like those 
in He-II or strongly type-II superconductors). 

Indeed, one can employ a vortex trial solution 



o*0 



R 



VR 2 



Z2 



VR 2 + r 2 



(43) 



only one com- 
ponent z\ = ■m = e 1 ^ survives on large distances. The 



which satisfies the boundary conditions 

e 



|A| 

solution (|4*3"| possesses a topological charge Q = 1, and 
it is just the skyrmion [l5|, [l6| for the N = 3 non-linear 
sigma model functional (i. c. the first term in the free 
energy fl2l ) 



-Eskyrr 



2l 

2 



-{d^n) 2 d 2 r. 



(44) 



Note, however, that here we use a different gauge. The 
solution (|43f corresponds to choice of the phases 



X = -¥>/2 = 0/2. 



(45) 



The parameter R is an arbitrary number: at any R the 
topological charge of the skyrmion is Q = 1 and its en- 
ergy 

-Eskyrm = ^Ps ■ 

However, the full gradient functional (|4"2")) contains the 
second term as well. This term leads to logarithmically 



E-> 



irp s log 



A 
R 



(46) 



where A is the minimal of the system size and (a very 
long) two-dimensional screening length X 2D = 2X 2 /d. In- 
deed, at r > R one finds At — r~ 2 , whereas at r <C R 
the polar angle 9 — > tt and, according to Eqs. (|4"ffl EES"]) 
the "vector potential" A^ is not singular anymore. It is 
evident from Eq. (gSJ) that one should choose R as large 
as possible in order to minimize vortex energy. The up- 
per limit is given by the anisotropy length L an defined in 
Eq. p6p . Thus the minimal energy of our trial solution 
can be estimated as 



Econt = 7rp s (log A/L an + C) . 



(47) 



where C ~ 1 . The energy of the continuous vortex (14"?| 
is lower than the energy of a usual singular vortex by a 
large amount 



F TP - 77 1 & 

losing — t^cont — ~^Ps *Og T^~\ 



(48) 



We emphasize, that the solution (|43|) does not provide 
the energy minimum but is just a trial function; the cor- 
rect non-singular vortex solution should have even lower 
energy and thus is more stable than the singular vor- 
tex. In the case of a continuous vortex the term Co in 
Eq. (|4"Tj) is zero and Q = <&/ <&0; whereas for a singular 
vortex Co = ±1 and Q = 0. 



C. Half-quantum vortices in the LOFF state 

The LOFF-like (or "stripe" ) state is realized at (3 a > 
and its degeneracy manifold is S 1 <£> S 1 . Indeed, the en- 
ergy minimum is realized when vector ft is parameterized 
as n — (cost^sin^jO), i. e. 8 = tv/2, thus there are 
two phase variables, <p and x> cf- parameterization ([571 
I38p . A usual singular vortex solution corresponds to ip = 
const and S\ = 2tt; according to Eq. (|40[) . in this case 
the "vector potential" = d^Xi an d it does not con- 
tain the second phase ip. Then a natural question arises, 
if some other vortex-like solutions are possible, due to 
the extended (with respect to the usual S 1 ) degener- 
acy manifold. Indeed, the same degeneracy manifold of 
the order parameter is realized in some of the p-wave 
superconductive states, leading to the existence of half- 
quantum vortice o 17 i 18 ' 19 . The reason for the existence of 
such an object is evident from the representation (|37[) : a 
sign-change of the order parameter due to the 7r-rotation 
of the phase x along some closed loop in real space can be 
compensated by the ±27r-rotation of the phase ip along 
(topologically) the same loop. 

We are not aware of any explicit solution for a half- 
vortex in the general case of an S 1 <E> S 1 degeneracy man- 
ifold. However, some progress can be made in the vicinity 
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of the symmetric point 5, where (3 a <C (3 S and the prob- 
lem simplifies a bit due to the presence of the "isotropic" 
spatial scales, f(T) < L < L an (T), cf. Eq. d3fij). In- 
deed, on such length-scales the problem can be treated 
within the gradient free energy, cf. Eq. (|35|) or Eq. (|42)l . 
The solutions with half-quantum of magnetic flux obey 
boundary conditions 6(00) = n/2 and Sooip — ±2tt, where 
via <5oo we denoted the phase increment along the large 
loop. Explicit form of these solutions may be (z~i = \[2zi) 



Zi = Jl- 



R 



Vr 2 + R 2 



R 



Vr 2 + R 2 
and 



R 



Vr 2 + R 2 



h = \ 1 



R 



y/r 2 + R 2 



e* 7 * (49) 



with an arbitrary parameter R. The variable 7 = ±1 
in the exponent of Eq. (j49f corresponds to the sign of 
the vorticity (the magnetic flux), whereas the first and 
the second lines in Eq. (|4"9"|) correspond to the cases of 
negative and positive projections 713 of the n vector in 
the center of the half-vortex. In terms of vector n these 
four solutions (|4"9"| are represented as follows: 



ni,2 




n 3 A 




An elementary calculation of the topological charge and 
of the magnetic flux associated with the vortex solution 
Eq. ((501) leads t0 



Q = - = 1 - 
$0 2 



(51) 



An additional binary variable which characterizes the 
half- vortex is the sign of the component n^(r = 0). 
Therefore, totally there are 4 types of half-quantum vor- 
tices. 

In the presence of a finite length L an the solutions (T4"9"l) 
make sense as intermediate asymptotics, as long as R <C 
L an , whereas at longer scales the anisotropic term cx (3 a 
modifies the solution considerably (a numerical solution 
would be necessary to determine the solution in this re- 
gion). Nevertheless, we can make an energy comparison 
between the singular vortex and a pair of half-quantum 
vortices even without an explicit solution including the 
anisotropic term. Consider two contributions to the free 
energy functional (|42p . The term with c^re, evaluated 
for the solution (|50[) . gives (the main contribution comes 
from large distances, r 3> R) 



Ex 



7T A 

4^ l0g i?' 



whereas the term with A^ contributes with about the 



same amount (since A^ is non-singular at small r < R) 
E A 



7T , A 

4^ log iT 



Both above estimates may contain subleading terms ~ p s 
which we do not control. Totally, the energy of a half- 
vortex is 



E1/2 = -^Ps 



log 



A 
R 



/2 



(52) 



where C1/2 ~ 1. The minimal energy of a half- vortex 
can be estimated from Eq. (|52[) with the substitution 
R ~ L an . Thus we find that the energy of two half- 
vortices coincides (up to the terms ~ p s which do not 
contain a large logarithm) with the energy of a contin- 
uous single-quantum vortex (I47p . found in the previous 
Subsection, and is certainly lower than the energy of the 
singular vortex by approximately the same amount as in 
Eq. (|48l) . This means that the half-quantum vortex is 
a fundamental topological defect of the stripe (LOFF) 
state, at least in the region relatively close to the sym- 
metric point S. 



V. PHASE DIAGRAM 

A. Stationary conditions for the helical phase 

Now we concentrate on the properties of the helical 
phase significantly below T C {H) 1 and determine the loca- 
tions of the phase transition lines £T, ST' and TO. This 
calculation is possible since |A(r)| 2 — const in the helical 
state, and thus explicit analytic equations determining 
A and the corresponding Q can be written without re- 
sorting to an expansion over small A. Evaluation of the 
thermodynamic potential in the helical state gives 



A 2 

n hel (A,H) = — - 

wv(eF)T 



(53) 



/■2-7T 



D 2 + A 2 



dip 
2^' 



where 



lu = lu + iH\ sin tp, 



(54) 



and H\ is determined in Eq. (|27[) . The equations de- 
termining A and Q are derived from the two stationary 
conditions 



dn hel = 2A 
OA U 

and 



„ 1 Jo 



A dip 



VJ> 2 + A 2 2tt 



(55) 



Oflhel _ VF 

dQ ~ T' 



■(e F )Tj2f(Hx,u;) = Q, (56) 



where we denoted f(H x , u) = -i v ^"^^- 

Reducing the integrals over ip in Eqs. (1551 EU) to com- 
plete elliptic integrals (see Appendix [T]) compacts the 



10 



stationary conditions to a two equations set 

1 9 I \ r F K(fc) 



J2 f(H x ,u) = 0, 



where the Jacobi modulus 



k = 



y/v* + (\H x \ + Ay 



and the function 
1 



K(fc) 



V^ 2 + (I^A| + A)2E(fc) 



(57) 



(58) 



(59) 



f(H x , = — [[^ + HI + A 2 ) _=_= 



(60) 



is defined through the Jacobi complete elliptic integrals 
of the first and second kind. At A = Eq. (|5T|) reduces 
to Eq. {281). 



The line £0 of stability of the BCS state (with respect 
to a formation of the helical wave) is determined by a 
simultaneous numerical solution of Eqs. (f57"l l58|) . taken 
at Q = 0, and Eq. (|6"4")) . This line is indeed a line of a 
second-order transition as long as the coefficient 6 > 0. 
Using Eqs. (I53[ l58|) . we compute the coordinates of the 
point T S CO where 6 = (cf. Eq. flM} for 6) as 
(if,T) = (1.547, 0.455)T cO . At lower temperatures 6 < 
and a first-order transition out of the homogeneous state 
takes place. Therefore TO is a boundary of a domain of 
the BCS state local stability. The actual first-order tran- 
sition line TO' between the BCS and some inhomoge- 
neous state consisting of many spatial harmonics lays at 
lower values of the magnetic field (Ho> < Ho = 1.76T c o). 

The homogeneous superconductor which exists on the 
left of the Lifshitz line is "gapless" for high enough 
temperatures. The spectrum of the particles is E p — 
■\/£ 2 + A 2 — XH sin (p p , therefore the minimum bound en- 
ergy of the Cooper pairs (the energy gap) turns to zero 
at H > A. The line of transition into the gapless super- 
conductivity H = A is marked in Fig. Q] with a dashed 
line. 



B. Lifshitz phase transition line CT 



The thermodynamic potential (|53[) in the helical state 
is a function of Q for any given pair (H, T) by virtue of 
Eqs. (H3 EH). At small Q Eq. dSH can be expanded in 
powers of Q as (terms of the order of a/vF have been 
neglected) 



fij,ei(Q) = n/iei(0) + aQ 2 + bQ i + cQ 6 , (61) 



where 

1 o 2 n 



2dQ 2 



A 2 sin <p dip 

( £2 + A2) 3/ 2 T' 



246 



d 4 E 



( o 3 n 
d A Vt \dZdQ 1 



(62) 



(63) 



and c > 0. The condition a = 0, b > determines 
the second-order Lifshitz phase transition line CT, which 
ends at the critical point T, where the coefficient 6 = 
changes sign, see Fig. [T] The condition a — may be 
simplified after reducing the integrals over ip in Eq. ([62]): 



UI>0 v 

where 



J = 



d 

2lu — J 

0u> 



A 2 



d 



dA 



J = 0, (64) 



C. Phase transition line ST' 

The second-order phase transition line ST' bounds the 
region of the helical phase from the high-i? side. Its 
position can be determined via the stability condition 
with respect to the additional modulations of the or- 
der parameter, of the form <5A(r) = 5v- q exp(-iqy) + 
8v q+2 Q exp(i(q + 2Q)y): 



SQg v = v^A(q)v, 
where v = (Sv- q , Sv* +2 q) and 

a( n \ _L ( G \p- G *-P+ F Xp _Fx- p+ 

[ '~U ^X F Xp- F x- P+ G Xp++Q G X - p _ +Q 



(65) 



(66) 



whith p± = p ± q/2. The Green's functions entering the 
matrix A are 

r 1 p " 2 ) -^ 2 + (c-s p ) 2 + a 2 ' 



Fx(u,p- 



Ae-^A 



2) ^ + (£ - Qp) 2 + A 2 



= -Fx ( -w,-p + | + Q 



(67) 



where a) is given by (|54[) . and for brevity we introduced 
a notation Q p = (g + Q) sm<p p /2. 

The matrix A has two eigenvalues £1(9) < £2(9), 
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where 

j_ _ 9-g + 9q+2Q 
U 2 



2n 



v / dlf . 1 fi- ^ 2_Q p 



.9-9 - gg+2g 



E 



w>0,A 



w>0,A 
2tt 



dly9 1 



o 4 Vci 2 + A 2 ti 2 + Q 2 + A 2 ' 

1 A 2 
T Vw 2 + A 2 a) 2 + Q 2 + A 2 ' 



(68) 



G\_ r ,_ n /')G\__ r ,_ 



w>0,p,A <J A,p-q/2 



p-q/2- 



^>o, p ,a Gx p+q /2 + QG A -p+q/2+Q and 



Here g- q 

9q+2Q = J2: 

/2 . The integrals (JH 
can be expressed through the elliptic integrals of the 
first and the second kind, as shown in Appendix l"2l 

The helical state metastability line ST' is defined as 
a set of points where one mode 5v becomes energeti- 
cally favorable: miriqei(g) = 0. Four equations are to 
be solved simultaneously in order to find the position of 
this line: two gap equations (|571 that determine the 
equilibrium A and Q, together with the two equations 
dq^i(q) = and e\{q) = 0, where €\{q) is determined 
in (f68|) . We call the new phase on the right of the ST' 
transition line a "three-exponential" state. 

Note, that the ST' line is an actual phase transition 
line out of the helical state if this transition is of the 
second order. Another possibility might be that a first- 
order transition occurs, which transforms the helical state 
into parity-even LOFF-type state, and occurs at slightly 
lower values of H at each T. Below we show that near 
T C (H) the phase transition is indeed of the second or- 
der. In order to demonstrate it, we evaluate terms of 
the eighth order in |Aq| in the Ginzburg-Landau func- 
tional. Parameterizing the order parameter spinor near 
the symmetric point as 



A_ 



Ae tx 



-itp/2 
Jip/2 



cos 9/2 
sin6>/2 



(69) 



the anisotropy part in the Ginzburg-Landau functional 
reads 



ants 
sn 



(«) =/3 a A 4 cos 2 



kA 8 cos 4 1 



(70) 



(note that the term cx |A| 6 is not allowed by symmetry). 
The sign of the coefficient k in front of the term cos 4 6 
determines the type of the transition near the T C (H) line. 
We demonstrate, by rather tedious calculations described 
in Appendix [31 that k > 0, which proves that the tran- 
sition is of the second order near the T C (H) line. 

Now we come to somewhat surprising situation. In- 
deed, according to the analysis of superconductive in- 
stability at the T C (H) line, on the right of the sym- 
metric point S the stripe ( LOFF-like) state is formed, 
which contains two harmonics with Q = ±|Q|. Such a 



phase preserves spatial inversion (in the plane) and time- 
reversal symmetry. Below the T C (H) line additional har- 
monics develop in such a phase, but they come in pairs 
±3Q, ±5Q, and still preserve spatial and time-reversal 
symmetry. On the other hand, the second-order phase 
transition line ST' separates the helical state (parity- 
odd) and a "three-exponential" state which has broken 
parity and time-reversal symmetry as well. It means that 
there must exist one more phase transition line, between 
the "three-exponential" phase and the parity-even stripe 
phase. This transition line should start at point S but 
will lie slightly to the right from the line ST', i. e. at 
higher values of the field. 

The fact that the points T and T' are different but 
close to each other is in favor of existence of a critical 
point K, on the line ST' (similar to the point T on the 
line CO) below which the phase transition from the heli- 
cal to the LOFF-type state becomes a weakly first order 
transition. 



VI. CURRENT AND ELECTROMAGNETIC 
RESPONSE IN THE HELICAL STATE 

A. The absence of the equilibrium current in the 
helical ground-state 

The oscillating space-dependence of the order param- 
eter like given by Eq. (fl6)) may lead to the hypothesis of 
a nonzero background electric current present in such a 
state. We will show here by general arguments that the 
electric current is in fact absent in the helical state: the 
condition of its vanishing is equivalent to the minimum 
of the free energy with respect to the variation of the 
structure's wavevector Q. 

The superconducting current can be written in the fol- 
lowing form: 

J s = £ T £ (^|±^G AX (^P + Q/2)- 

u?,p,A 



= A ,-p+Q/2 

dp 



G AA (-^,- P + Q/2) 



where 



G A a(^,p + Q/2) 



*W + eA,-p+Q/2 

" cl> 2 + £ 2 + A 2 



(71) 



(72) 



is the electron Green's function in the helical state; e AjP 
is the spectrum of the free electron. For brevity we used 
notations 



LJ = LU + i 



eA,p+Q/2 - £A,-p+Q/2 



$A 



eA,p+Q/2 + e A : -p+Q/2 
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The thermodynamic potential in the helical state 



T 



hel 



log [o) 2 + ^ + A 2 ] 



(73) 



is evaluated explicitly for an arbitrary spectrum e\, p of 
the electron, due to |A(r)| =const. 
The derivative is 



dQ 
4 ^ 



A,p + Q/2 

9p 



9fx 



-p+Q/2 



(74) 

(-*£>+&) 



■3 



A 2 



Comparing this stationary condition (|74[) with the ex- 
pression for the current (f7"Tj) . we notice 



2e gOfee, 



(75) 



Thus, a direct calculation of the superconducting current 
j s shows that in any order of cx/vf the current in equi- 
librium is zero. This result apparently contradicts the 
statement made by Yip2£, who studied basically the same 
model as the present one, and found a non-zero super- 
current in the presence of a parallel magnetic field. The 
resolution of this paradox will be presented in Sec. IVIII 
below. 



B. Electromagnetic response in the helical state 

We calculated the electromagnetic response function 
Sja/SAp = ——n" 13 for the helical state using the stan- 
dard diagram methods: 



Jx 



''-TJ2 (G 2 +F 2 )co S 2 <p(^~\ a y A x , 



u>.p,X 



^ = 7 T E (G 2 + F 2 ) sinV - Aa) A y , 

ui.p,X 

(76) 

where G = G(w,p + Q/2) and F = F(w,p + Q/2) are 
correspondingly the normal and the anomalous Green's 
functions in the helical state. We found that 



d 2 n 



h dQ 2 



(77) 



i. e. proportional to the parameter a from the expan- 
sion of the thermodynamic potential (foTj) . Thus on the 
Lifshitz line CT there is no linear supercurrent in the 
direction perpendicular to the magnetic held. The com- 
ponent n xx does not vanish anywhere in the helical state 
region and is of the order of n s of the BCS state: 



4m 

IF 



(f 




f i dn hel \ 


d 2 n hel \ 




dA 


V A OA J 


dQ 2 J 



(78) 



This anisotropic behavior of the superfluid density tensor 
is in contrast with the one found in the classical LOFF 
problem, where n s was shown to vanish in the whole he- 
lical state; the difference is probably due to the fact that 
in our problem the direction of Q is fixed by the applied 
field h, while in the case of a ferromagnetic supercon- 
ductor it is arbitrary. The obtained behavior of the n"" 
tensor indicates a strongly anisotropic electromagnetic 
response of the surface superconductor near the Lifshitz 
line CT. 



VII. "WEAKLY HELICAL" PHASE AT LOW 
MAGNETIC FIELDS 

A. Transformation of the uniform BCS state into 
long-wavelength helical state 

The thermodynamic potential of the helical state in the 
form of Eq. (|53|) was obtained while neglecting the term 
aQ small in comparison with the term vpQ in the en- 
ergy of the electron. Within this approximation the ther- 
modynamic potential was symmetric under the change 
Q — > —Q and hence the expansion (|61|) contained only 
even powers of Q. The account for the first-order term 
in a/vp in expressions of the type of Eq. ([53"j) leads (see 
below) to an appearance of the term r\Q in the thermo- 
dynamic potential (foTj) . This term results in the trans- 
formation of the homogeneous BCS state (situated on 
the left of the CT line) into a weakly helical phase with 
a small wave-vector \Qhei\ = ^r-, first found in fnl for 

V F 

small values of H. Mathematically, the terms of the order 
of a/vp are considered via taking into account the differ- 
ence in the density of states in Eq. |53|) : v{ep) — > V\{ep). 
Then the second stationary condition ([58]) is changed no- 
ticeably: 



dQ 



K*F)T£(Y + A f)/( ff A,w)=0, (79) 



to , A 



where f(H x ,uj) is defined in Eq.®. Now the Q = 
solution never (at any H ^ 0) provides a minimum for 
the superconducting energy: 



hi I 



dQ 



3=0 



av(ep)T^2f(H,cj)^0, (80) 



and the expansion of the thermodynamic potential in 
powers of Q contains the linear term: 



(81) 



which clearly indicates that the equilibrium wave vector 
modulating the order parameter on the left of the CT 
line is always non-zero, although small as a/vF <C 1: 



Q 



hel 



rj_ 
2a 



2aH 



(82) 




FIG. 7: The dependence of the pairing wave vector Q upon 
the magnetic field for the two values: a/vF = 0.1,0.01 (and 
a particular temperature T = 0.54T c o). The curves shown 
by bold and thin lines correspond to neglecting and taking 
into account in the SC energy of the helical state the small 
terms o/vf- The dashed lines correspond to the approxima- 
tion (JH2J. 



y 


\ {uj2+H * +A2) K(k) VE(fc)" 









where k is the Jacobi modulus (|59|) and V = 
Y^ 2 + (H + A) 2 . In the limit H -> Eq. d52J) can be 
expanded, resulting in the dependence linear in 7?: 

Qfcei = -2aH/v%. (83) 

The formula (|83|) is applicable for small enough magnetic 
fields; in particular, for T w 0.5T c o the field range is 
limited by H < 0.5T c0 , cf. Fig. At higher fields the 
dependence of the pairing wave vector upon the magnetic 
field becomes non-linear but still can be approximated 
by Eq. fl52J), until H < 1.3T c0 , cf. Fig. Within the 
above range of magnetic fields the expansion ([81]) of the 
thermodynamic potential is still applicable. 

When the magnetic field is further increased, the sys- 
tem approaches the transition into a short-wavelength 
helical state discussed in Secs.V and VI above. The ma- 
jor role of the small term r]Q is to broaden the CT line 
of the second-order transition from the uniform to the 
helical state into a narrow crossover region. Now the de- 
pendence of the pairing wave vector upon the magnetic 
field should be obtained numerically by solving the sys- 
tem of the two self-consistency equations (f5*o*j) and (f5*8")) . 
while keeping also the term oc a/vp via the substitution 
of v{ep) by i>{ep)(l + Xa/vp). The corresponding results 
for Q(H) are shown by full thin lines in Fig. for two 
values of a/vp. 
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FIG. 8: A superconducting film wrapped in a cylinder (cyclic 
boundary conditions) . This geometrical configuration enables 
a current to flow in equilibrium in the ground state of the 
helical phase when a parallel to the axis magnetic field is 
applied. 

B. Cancellation of the ground-state current and a 
spin-orbital analog of Little-Parks oscillations 

The gradient of the phase of the condensate wave func- 
tion determines the density of the superconducting cur- 
rent 

= ^n s Q hel , (84) 
2m 

where n s is the density of the number of superconducting 
electrons, e = — |e| is the charge of the electron, and m 
is the electron true mass. The expansion of Eq. ([77]) for 

2 

weak magnetic fields gives n v s v — , -^ IL iy(ep)(l — Y(T, A)). 
In this limit the density of the current in y-direction, 
induced by the superconducting phase gradient, reads 

jW = -ev(e F )a(l - Y(T, A))H. (85) 

The supercurrent given by Eq. (|85[) is not the only con- 
tribution to be considered. In fact, Yip2£ have considered 
the same model and found that a weak current propor- 
tional and perpendicular to the magnetic field flows in 
the uniform BCS state. Indeed, a calculation of the cur- 
rent in the BCS state gives a non-zero value (coinciding 
with that in [H) 

4 2) = T E = ev{e F )a{l - Y(T,A))H. 

LJ,p,A 

This second contribution leads to the presence (due to 
the Rashba term in the Hamiltonian) of the anomalous 
contribution to the electric current. 

In the true ground-state with the weak helical modu- 
lations, both contributions to the current, Eqs. (|55l 

sum up to produce a perfect zero: jlp + jy 2 ^ = 0, as 
they should do according to the general proof given in 
Sec. I VI Al above. Thus we see the resolution of the con- 
troversy with the result (2(|: a uniform BCS state con- 
sidered by Yip does produce a supercurrent under the 
action of a parallel magnetic field, but this state is not 
the ground-state. Instead, the ground-state is realized as 
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in-plane magnetic field 



FIG. 9: The sawtooth dependence of the superconducting 
current j, flowing around a cylinder, upon the magnetic field 
H. The amplitude is j ma x = ^n s /R, where R is the cylinder 
radius. The period is H — Tiv%/(2aR). The strictly linear 
dependence, shown on the figure, takes place in the case H <C 
T c0 , T = only. 



a weakly helical state with a zero current. In fact, when 
a parallel magnetic field is applied, the current does not 
flow, but a phase difference A% = LQhei is induced at 
the edges of the superconducting film in transverse to the 
field direction. 

However, interesting "traces" of the spin-orbit-induced 
supercurrent could possibly be seen, if the superconduc- 
tor is wrapped in a cylinder and a magnetic field is ap- 
plied along the axis. Then a current may flow around the 
cylinder (see Fig. [SJ, due to the quantization condition 
8\ = 27m, n = 0, ±1, ±2, which does not allow an ar- 
bitrary phase shift along the closed loop which encircles 
the cylinder. The total current will be given by a sum of 

(2) 

the "spin-orbit current" y s and the current due to the 
gradient of the quantized phase: 



Jquant 



eh 
— ? 
2m 



Ihel 



2im\ 
2^RJ 



where R is the cylinder radius and Qhei = —2ctH/v F . 
The integer number n is determined as a function of the 
magnetic field in a way to minimize the current: n = in- 

2 

teger part oi[H/H ], where H Q = The current will 

vanish at the field values H = nH a , n = 0, ±1,±2, ... 
only. The maximal value of the current is equal (at 
T = 0) to j m ax — 7r*-n a /R. The dependence of the 
current on the magnetic field is of the sawtooth form, 
as shown in Fig. [9j in the ideal case of zero tempera- 
ture and absence of impurities. The predicted oscilla- 
tions of the current are, on first sight, similar to the well- 
known oscillations in a superconducting cylinder (with 
the period Blp = Tic/eR 2 in terms of a real magnetic 
induction), which goes back to the early stages of the 
superconductivity studies. 21 However, the correspond- 



ing oscillation periods differ by orders of magnitude: 
H /H LP ~ {vf/ol) ■ (k F R) > 1. We note also, that 
in a real experiment it should not be necessary the su- 
perconducting film to be wrapped in a cylinder; it would 
be sufficient to fabricate a heterostructure where a thin 
film with a SO coupling would serve as a weak link in- 
troduced into a SQUID-type loop, which would allow to 
control the superconducting phase difference between the 
film edges. 



VIII. PHASE DIAGRAM IN THE PRESENCE 
OF NON-MAGNETIC IMPURITIES 

In this Section we study the effects of potential impuri- 
ties upon the superconductive instability in the presence 
of the Rashba term and the parallel magnetic field. We 
consider a standard model of short-range weak impuri- 
ties with a potential u(r) = uS(r) and density ni mp ; they 
are characterized by an elastic scattering time r, where 
r -1 = 2irni m pU 2 v(£p). It is assumed that the impurity 
scattering rate is weak with respect to the Rashba split- 
ting, 1/t -C apF, whereas it can be both weak or strong 
in comparison with T c . 

The interaction between the electrons and the impurity 
is described by the following Hamiltonian written in the 
chiral representation: 



Hi, 



V 



ue 



-»(p- 



(87) 



where the matrix 



M A/1 (p,p') = i(l + A Me i ^-M) 



appears due to the transformation from the spin to chiral 
basis. The semiclassical electron Green's function of the 
chiral metal in the presence of non-magnetic impurities 
is 



G? A (w,p) = 



1 



£a(p) - XHsinipp + ^sgnw 



(89) 



where £a(p) is determined in (|22"|) . 

The electron-electron vertex in the Cooper channel is 
then given by non-crossing diagrams shown in Fig. 1101 
where G denotes the Green's function (f8"9")) . Note that 
the chiralities of both electrons which belong to the same 
"Cooper block" coincide (otherwise, a diagram would be 
smaller by a factor (apFr)^ 1 <C 1). 

An analytical expression for the impurity line is (cf. 
Eq. (|55|) for the matrix elements) 



V\n((p p ,ipp>) 



1 



8TTTu(eF 



■(1 + A/ze^-^p') 2 = (9°) 



Aktv^f) 



p' ( A// + sin ip p sin ip p > + cos ip p cos (p p > ) , 
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FIG. 10: The electron-electron vertex in the Cooper channel 
in the presence of non-magnetic impurities. 



where A and \i are the chiralities of the left and the right 
blocks around the impurity line. In fact, the last term 
in the r.h.s. of Eq. ([50)1 can be safely omitted since its 
contribution vanishes after integration over the momenta 
in the product of and the Cooper block, Eq. (|9Tj) be- 
low. Integrating over £\ a block of two Green's functions 
(a "Cooper block") in the Cooper channel gives 

C\(uj, sin tp p ) = 

= v x (e F ) J d£ G x (v, P + |) G A (-w, -p + ~ 



iu: — H\ sin ip p 



(91) 



where u> — lu + sgnw/2r, H\ is determined in (|2T[) and 
v\(^f) = v(e F ){l + Xa/vp). 

The superconductive transition temperature is deter- 
mined by the usual condition ?7(0) _1 = C, where C is the 
sum of all ladder diagrams with n — 0, 1, 2, ... impurity 
lines (see Fig. ITU]) : 



w>0n=0 A„ 



2tt 



2tt 



xC A >,sin^ p JA„e^", (92) 

where L A ' (u>, ip Pn ) is the expression for a ladder diagram 
which includes a left vertex L° Xq {to, tp po ) = Aoe~ I¥,p o, n 
"Cooper blocks" and n impurity lines; the factor A„e lVp ™ 
corresponds to the right vertex of the diagram in Fig.fTOl 
In order to sum up the whole Cooper ladder, it is use- 
ful to employ a recurrent relation between the ladder di- 
agrams of the n-th and the n + 1-t order: 



^: + >^p n+1 )=E 

A„ 



d(f Pn 
2tt 



L \ n (u,ip p J x 



x Ca„ (w, sin ip Pn )V\ n \ n+1 {ip Pn , <p Pn+1 ). (93) 

The form of Eq. (j9"0)) helps to identify an Anzats 

L n Xn (w, cp Pn ) = {l° n (An, w) + J n (A„, w) sin^pje"^- (94) 

for the solution which is consistent with the recurrent 
relation (|9"3"|) . After substituting the Anzats (|94p in 
Eq. (|9"3")l , we encounter integrals 



f.27T 



diy9 i sin 3 tp 



4r J q 2n iuj — H\ sin ip 



, where j = 0,1, 2.(95) 



Then equation (|93[) can be rewritten in the matrix form 

i n+1 = m ni (96) 

where we define a 4- vector 

£ = (#(+),£(+)>#(-),£(-)) (97) 

and a 4 x 4 matrix 



R = 



4 /| /i 
ii 7i /i p 



(98) 



with the three "block integrals", Eq. (|95|) . calculated as 

1 1 




(99) 



Now we can proceed with the calculation of the Cooper 
ladder. We substitute the Anzats and the Cooper 
block (|9~lj) in the Cooper ladder ([92]) . and obtain 



C = 4^( ej ,)r^^^A^(A)/ A + ^(A)/i) = 

aj>0n=0 A 

= 47TTiy(e F )r^/ T • (l-i?) _1 -fo, (100) 



w>0 



where we used the definition of the integrals (|95[) and the 
4- vectors (|97p. and we summed up the geometric pro- 
gression X^^Lo^™ = (1 — -R) -1 - We also introduced 
a vector = (4, J]_, — — J|) and used a notation 
^ = (1, 0, —1, 0), corresponding to the definition of 
L\ (o>, ip Po ). At this point we have as well neglected the 
difference v+ ^ V—. 

Evaluating the scalar product in Eq. (|100[) . we find the 
Cooper ladder C and the equation for T c : 



1 



v{e F )U{0) 



UJ>0 



(101) 



where the kernel is 

(J° -Ilf 



K{u) = 4r- 



(l-(J°+J°)) [l-(/2 +/ 2)]_ (/ l _J1)2- 

(102) 

Equation (|101[) for T c (/i) was solved numerically; the 
phase transition lines are shown in Fig. [11] for a number 
of impurity scattering strengths l/2rT c o, interpolating 
from a clean to dirty limit. It is seen that in the clean 
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Phase transition lines for different strengths of impurity scattering: 



0.0, 0.5, 1.2, 2.1, 3.2, 4.5, 6.0, 7.7, 9.6, 11.7. Lifshitz points are shown by crosses. 



1/2tT c0 



limit T c qt 3> 1 the impurity scattering decreases the crit- 
ical parallel magnetic field (in the clean limit it is given by 
H p v = A(0) first found in and simultaneously 

pushes the position of the C point to higher values of H 
and lower values of T. As a result, both short- wavelength 
inhomogeneous states disappear from the phase diagram 
at t -1 > 9T c q. Note a large numerical factor 9 which 
enters the definition of the "dirty" regime in the present 
problem. In the "dirty limit" 1/r > 10T c o the kernel 
K(uj) simplifies to 



dirty limit the kernel reads 
K(lo) = - 



2H 2 t + v 2 f Q 2 t/4 + 2raQH ' 



(105) 



Maximization of the Cooper loop with the kernel (|105[) 
with respect to the pairing vector Q leads to a non-zero 
momentum of the Cooper pair, equal to 



!hel 



(106) 



2Hh 



,Q 2 t/A 



(103) 



The kernel (jl03[) is maximal at Q = 0, i. e. in the 
dirty limit one obtains a homogeneous superconducting 
phase (see below, however). The zero-temperature limit 
of the paramagnetic critical magnetic field H p (0) can be 
easily obtained with the use of Eq. (|103|) as 



for all values of the magnetic field. Thus the terms a/vp 
transform the homogeneous superconducting state into a 
weakly inhomogeneous helical state, analogously to the 
clean case. Note, that the small wave vector modulating 
the order parameter in the "dirty" limit, Eq. (|106[) . is 
twice larger than it is in the clean case for weak magnetic 
fields, Eq. (|83[). Note, however, that the result (|106|) was 
obtained near the transition line T C (H) only. 



H p = 



4reT ' 



(104) 



Thus, in the dirty limit 1/r 3> T c o the paramagnetic 
critical field grows with increase of disorder. 

The above results in this Section were obtained in the 
main order over the parameter a/vF- The linear in a/vF 
terms can be included in the same way it was done in 
the preceding Section, which lead to the substitution of 
Xcc/vf) in the kernel (|102p . Then in the 



3(1 



IX. FLUCTUATIONAL EFFECTS NEAR THE 
T C (H) TRANSITION LINE 

The corrections to the mean-field approximation, em- 
ployed in this paper, are usually of the order of T c /ep 
for a clean 2D superconductor: the actual transition is 
of the Berezinsky-Kosterlitz-Thouless vortex depairing 
type, and the transition temperature is shifted down- 
wards by a relative amount T c /ep < 1. In our system 
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the fluctuations are enhanced strongly around the special 
points C and S, where an additional analysis is needed. 

We start from the C point and recall that near the 
whole CT line the component n v s v of the superfluid den- 
sity is suppressed, cf. Eq. (177|) . proportionally to the 
coefficient a from Eq. (fr?2")) . Within the approximation 
used in Sec. IVII and above, this reduction factor could 
be arbitrary small, since n v s v vanishes exactly at the 
CT line. However, as it was explained in Sec. IVII Al 
an account of the sub- leading terms cx a/vF trans- 
forms the CT line of the Lifshitz phase transition into a 
crossover region. The minimal value of the second deriva- 
tive d 2 n he i(Q)/dQ 2 and thus of the ratio nf /nf x then 
scales as rj 2 ^ 3 cx (a/v^ ) 2 / 3 . In a superconductor with an 
anisotropic tensor of the superfluid density, the effective 
"rigidity modulus" is controlled by the geometric average 



The strength of the phase fluctuations is thus 



larger near the CT line by a factor (vp /a) 1 / 3 . Therefore, 
the fluctuational reduction of the transition temperature 
near the C point is enhanced by the same relative factor 
{vp/a) 1 / 3 > 1, cf. Fig.HU 

Another mechanism of a fluctuation enhancement is 
effective near the S point due to the extended U (2) sym- 
metry of the order parameter. Exactly at the symmetric 
point S the order parameter spinor spans the sphere S 3 
and is equivalent to a four component unit vector, see 
Eq. ([33)) . The thermal fluctuations of the classical O(N) 
nonlinear vector model in 2D space have been studied by 
Polyakov [22l| . On a large length scales L £(T), where 
£(T) is the temperature-dependent superconductive cor- 
relation length, the evolution of the dimensionless cou- 
pling constant g = T {i^fiA is governed (for N = 4) 
by the renormalization-group equation 



dg 
dX 



-9 



(107) 



where X = log(L/£(T)). Deviations from the symmet- 
ric point are measured (cf. Sec. IVI|) by the anisotropy 
parameter (3 a which enters into the free energy in the 
combination 



(3 a (T, H)\A\\M 2 + N 2 - N 2 - Ml) 



■2\2 



(108) 



It is easy to show the following—: that the anisotropy 
parameter (3 a satisfies a renormalization group equation: 



dA _ _ig 

iv — Pa ■ 

dX TT 



(109) 



Running solution of Eqs. (|107l 1109"]) can be easily found: 
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FIG. 12: The phase diagram with the physical Tbkt(H) tran- 
sition line shown by a dotted line, for the values T c o /ef = 0.02 
and a/vF = 0.34. The effect of the enhanced fluctuations in 
the regions near the points L and S is seen as a shift of the 
Tbkt(H) line in the direction of lower temperature and field, 
with respect to its mean-field position. The mean-field phase 
transition line T C {H) is shown by a solid line. 



(a renormalization of the coefficient (3 S is absent within 
the same approximation). The infra-red cutoff length 
L max (T) for the renormalization flow solution (|110l llllj) 
coincides with the length L an defined in Eq. ([36]) . Finally, 
we obtain the renormalized parameters 



2Trn /# 

— r^iog 
irn 



Pa 



and 



&a = Pi 1 



2Tm 

irh 2 n° g {f3 a 



(112) 



(113) 



The above formulae are valid as long as the effective 
renormalization group "charge" g{X raax ) is small com- 
pared to unity; within this domain a strong suppression 
of (3 a is still possible. Qualitatively, the result of the 
C/(2)-fluctuations is twofold: i) the "nearly-isotropic" be- 
havior extends to a wider region around the S point, and 
ii) the transition temperature occurs to be additionally 
suppressed within the same region. 

To conclude this Section, fluctuations lead to the defor- 
mation of the phase transition line T C (H) in the vicinities 
of the C and S points, as shown in Fig. [12J 



X. CONCLUSIONS 

In this paper we found the phase diagram of a surface 
superconductor with a relatively strong Rashba splitting 
of the electron spectrum: T c <g; app <C €f, in the pres- 
ence of a parallel magnetic field, both in clean and disor- 
dered cases. In the clean limit, we demonstrated that in 
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the lowest approximation over the spin-orbital parameter 
a/vp the phase diagram is universal and contains the un- 
usual "helical" phase with the order-parameter modula- 
tion oc exp(iQr), as well as a usual BCS-type phase and 
a LOFF phase with a sinusoidal modulation at higher 
magnetic fields. The absolute value of the modulation 
wavevector is typically Q ~ g^Bh/hvF- Non-magnetic 
impurities tend to diminish the region where these mod- 
ulated phases exist, and eliminate them completely in the 
dirty limit. Once sub- leading terms of the order of a/vp 
(which are due to the different chiral subband DoS) are 
taken into account, the uniform BCS state becomes un- 
stable and transforms into a "weakly modulated" helical 
phase with Q s iow ~ agfish/hvp. A weakly modulated 
phase is stable with respect to disorder, as its origin can 
be found in the symmetry of the problem^ 

In the clean limit, we were able to determine the posi- 
tions of the phase transition lines between the BCS, the 
helical and the LOFF states on the (H, T) phase dia- 
gram. It was found possible to implement this calcula- 
tion mainly analytically, due to the fact that the absolute 
value of the order parameter is constant in the helical 
state. We have shown that the thermal phase fluctua- 
tions are enhanced strongly around both these transition 
lines, leading to local "caves" of the T c (h) transition line 
near the end-points C and S. We expect that a fluctu- 
ational correction to the conductivity of the Aslamazov- 
Larkin type 2 ^ should be strongly anisotropic near the £ 
point, which could be one of the signatures of the pro- 
posed phase diagram. 

We have shown that an electric current does not flow in 
the equilibrium state of our system at any parallel mag- 
netic field (as long as the system is infinite), contrary 
to the statement made in [20(. However, if the system is 
considered in the finite-stripe geometry with the periodic 
boundary conditions, an oscillating (as a function of the 
parallel magnetic field) current is expected, which is an- 
other special feature of a superconductor with a strong 
Rashba coupling. 

A new type of a vortex-like topological defect was pre- 
dicted to exist in the region near the special S point of 
the phase diagram, due to the extended U(2) order pa- 
rameter symmetry realized at this point. Contrary to 
the usual singular vortex, this new topological defect is 
non-singular in the sense that the amplitude of the or- 
der parameter does not vanish in the vortex core. We 
showed that the energy of the non-singular vortex is lower 
than that of the singular vortex in some finite region 
around the S point. The nature of an elementary vortex 
defect differs between the helical and the stripe phases 
(in the regions to the left and to the right from the S 
point): whereas in the helical state a vortex carries an 
integer flux <!> = n$o and an integer topological charge 
Q = n, these quantum numbers are half-integer in the 
stripe state. Physically, these half-vortices correspond 
to the presence of a vortex-like defect within only one 
(either A + or A_) component of the order parameter. 

Note added. — When this paper was nearly completed, 



we learned about a recent preprint by Agterberg and 
Kaur— , where the density of states effects (due to v + ^ 
V- ) upon the phase diagram of the Rashba superconduc- 
tor were studied in the clean limit by means of a numer- 
ical solution of the Eilenberger equations. The main em- 
phasis was made on the 3-dimensional systems there, thus 
the comparison is not straightforward. Qualitatively, the 
results they obtained seem to be in agreement with ours. 
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APPENDIX 



1. Transformation to elliptic integrals 



Consider the integral 



h 



dip 



{uj + iHx sin ipf + A 2 



which enters Eq. (jB"5")) , and note that f Q 2 * f (simp) dip 



(114) 



4 Jq f(cos2ip)dip for any function /. Using then the 
identity cos lip = cos 2 ip — sin 2 ip and the substitution 
ip = arctani, we transform Eq. (|114|) into 



h=2 



dt 



Jo .jAe^t A + 2(uj 2 + H 2 + A 2 )t 2 + Aer^ ' 

(115) 

where we denoted Ae 1 ^ = (u> — iH\) 2 + A 2 , and 

A = ^(u 2 + Pa| + A) 2 )(iv 2 + Pa] - A) 2 ). (116) 



The integral (|115[) is rapidly convergent, thus it is possi- 
ble to make a transformation of variables in the complex 
plane r = te l ~ and to reduce Eq. (|115[) to 



h =2 



dr 



o VA^/l + (a + a- 1 )T 2 + t 4 



where 



V + (|g A | + A) 2 ) 
(oj 2 + (\H\\ - A) 2 ) ' 



(117) 



(118) 
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(in the course of transformation from Eq. (| 1 1 5|) to 
Eq. (|117[) we rotated the path of integration over r by the 
angle ip/2). Making an inverse substitution r = tan<^ in 
Eq. (TTT7)) gives 



2 r 12 

Ah 



dip 



(119) 



cos 2 2ip 



Following the definition of an elliptic integral of the 

first kind K(fc) = J£^(l ~ k 2 sin 2 tp)~ 1 / 2 dtp, we rewrite 
Eq. (fTT9|) as 



h 




(120) 



The last expression in (|120p leads directly to 
the gap equation (|57[) . with A and a given by 
Eqs. (fTTo)) and 1(115)1 . 



2. Equation for the ST' line 

Here we present the transformation of Eq. (|68|) , which 
determines the helical state metastability line ST' , to the 
complete elliptic integrals of the first and the third kind 
(K and II): 



h = 

and 

h = 



dtp- 



1 -1 + (uj + iX sirup) 2 



o 4Vw 2 + A 2 Cj 2 + A 2 + {X sin <p) 2 
zK(k) ~ Zl n{l u k) - z 2 U(l 2 ,k) 



(121) 



2tt 



1 



A 2 



o 4V^ 2 + A 2 Cj 2 + A 2 + (X sin <p>) 2 
yK(k)-y 1 n(h,k)-y 2 n(l 2 ,k) 



(122) 



where Q — uj + iH sintp. In Eqs. (rf2T)[T22)) we introduced 
notations 

(A-g + »c)(A 2 -(g + ^) 2 ) + A 2 ) 
U ' {A + H + ioj){A 2 + H 2 + uj 2 - X 2 - 2i$y 

Z 2 (S) = Ji(-S); 



-2- 



Zl (S) 



AH(2iuX + A(H + 2X)) 
(A + iw) 2 A 2 ' 
(A-iT + iw) L(S) 



2 (A + iu f X 2 (A 2 -{H- iu f + X 2 ) S 



*i(-S); 

A 2 ff 2 
(A + iw) 2 X 2 



w(s) 



iA 2 (A-H + iuj) M(S) 



2 (A + iw) 2 X 2 (A 2 - (iT - iw) 2 + X 2 ) S ' 
Ite(S) = yi(-S); (123) 



where 



L(S) = -iA 4 H(H + X) 2 - 
2u 2 X 2 (H + X)(luX + S) + 
2iAuX(H + X)(H -iuj + X)(luX + S) + 
A 3 (# + X) 2 (-iH 2 -Hlu + 2ujX + iX 2 + S) + 
A 2 (ff + A)(H 2 - iHuj + HX + iujX)(2ujX + S) 
uX 2 A 2 {H + X) 2 , 



S = yJ-A 2 B 2 + (A 2 +lu 2 )X 2 , X = 



q + Q 



(124) 



M(S) = AH 3 {A + H -iu)- 

- HX 2 (A(A + 2iT) - i(A - i3> + 2w 2 ) + 
+ {A + iuj)X 4 + 

+ (H 2 (i{A + H)+lj) + (-i(A + H)+u;)X 2 )S. 

(125) 



3. The eighth order terms in the Ginzburg-Landau 
functional 

Here we present a calculation of the coefficient n which 
is in front of the eighth-order anisotropic term in the free 
energy, Eq. ((70)l . As it was mentioned, the sign of k 
determines the type of the helical state instability. 

The terms of the eighth order in A in the Ginzburg- 
Landau functional arc 



F. 



(8) 



-A s 



- 8 m 



1 A 4 



(D 1 -D 2 + D 3 /2) 



(M 2 -M 2 ) 4 + £ 3Q , (126) 



where all but the last term originate from the expan- 
sion of the free energy in powers of the basic ampli- 
tudes u and v, and notations a = D\ + D 2 + -D3/2 
and b = 6D1 — D 3 are introduced. The last term E 3 q 
arises due to the nonlincarity-induced additional harmon- 
ics with momenta ±3Q and amplitudes u± 3 q and v± 3 q; 
explicit form of E 3 q will be specified shortly. Actually 
we are interested in terms which contribute to the coef- 
ficient k in front of the cos 4 9 term in Eq. ((70)1 , thus the 
terms in the first line of Eq. (|126|) are irrelevant. 

The coefficients f 1,2,3 come from the 8-vertices dia- 
grams shown in Fig. [T3J The corresponding analytical 
expressions are summarized below (we introduce nota- 
tions A=T£ WiP ,aA): 
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Di 
D 3 



1 s-iA f-ii 

2 P+T -P+T -P—T 

^ ,G 2 



/"to f 'il f < 

^ p+ Q^- p+ qG_ p _ 3 _qG p+5 _q, 



G, q G_ 



3(1 



G. 



+ 



S~lZ f~i£ ri s < 

G p +# tT -p+f G -p-^ G p-^ G p+^' 



(127) 



where G are the normal metal Green's functions. 
The integration is first done over e?£ = d(p 2 /2m — 
Ep), then over dip by means of a generating function 
1/^/uj 2 + {XH ± nQ) 2 , where n = 1,3,5. The obtained 
analytical expressions are evaluated for the values T$ = 
1.779T c0 , H s = 0.525T c0 and Q s = 2.647T c0 at the S 
point of the T C (H) line. The summation over uo is per- 
formed numerically, resulting in 



where 



D 1 = 0.00106053, 



D 2 = 0.00152674, 
D 3 = -0.00120067. (128) 



Now we turn to the evaluation of the "induced" term 
E 3 q in the free energy (|126[) . It appears due to a genera- 
tion of the third harmonics e Q in the order parameter. 
The amplitudes of the third harmonics u 3 q and v 3 q are 
small and proportional to the third power of the basic 
amplitudes u and v. The "induced " term reads 



E 3Q = V + AV + {X+V + V + X), 



(129) 



where the matrix A and the vector X contain numerical 
coefficients Z,-, which are determined by means of evalua- 
tion of the loop diagrams with two, four and six vertices 
(shown in Fig. [T4|) : 



A 



Zi + l 3 \u\ 2 + l 4 \v\ 2 l 5 uv 

i 5 u*v* z 1 + z 4 M 2 + z 3 M 2 



V + = (u* 3Q ,v. 3Q ), X+ = (x*,y), 

x = l 2 u 2 v* + l 6 u 2 v*\u\ 2 + (Z 7 + l$)u 2 v*\v\ 2 ; 

y = l 2 v 2 u* + l 6 v 2 u* \v\ 2 + (l 7 + h)v 2 u* \u\ 2 . (130) 

The fifth and higher harmonics, as well as the higher pow- 
ers of the third harmonics, do not contribute to the eighth 
order in A. The amplitudes u 3 q and V- 3 q in Eq. (|129p 
are found from the minimization of the superconducting 
energy: d ^ Q = 0, which gives 

E 3 q = -X+A^X = 

d 2 -(\u\ 6 \v\ 2 + \u\ 2 \v\ 6 )+d 3 -\u\ 4 \v\\ (131) 



d 2 = 



l 2 {-2hh+hh) 



Z 2 



2Z 2 (-2(Z 7 + Z 8 )Z 1 + Z 2 (Z 4 + Z 5 )) 



(132) 

The analytic expressions for the diagrams shown in 
Fig. [14] are (we introduce notations U = T A Zi) 

il = ^^E G p+ sg_G_ p+ 3_g_, 
w,p,A 

= 7>G p+ aG_ p+ QG_ p _3QG p+ sQ, 

1 3 = G_ p+ qG 2 p _3qG_ p ^3Q, 

h = G_ p+ SQG 2 p _3QG_ p _3Q, 

h = G_ p+ QG p _3QG_ p _3Q_G p+ 5Q, 

<6 - G p+ QG p+ ^G p _^G^ p _3Q, 

h - 2 G P+ § G - P+ %G p _^G_ p _ 3 _z, 

k = G p+ ^G p+ qG_ p+ qG p _3qG 2 _ p zq. 

(133) 

The calculation of the above integrals in the symmetric 
point S leads to numerical values 

Zi = 0.674254, Z 2 = 0.103447, 
h = 0.191039, Z 4 = 0.170477, 

h = 0.206894, Z 6 = 0.0325726, 
l 7 = 0.00491175, l s = 0.021198, (134) 

which we substitute in Eq. (|132p . Finally for the coeffi- 
cient k we obtain a positive value: 

k = D 1 - D 2 + As/2 - d 2 + d 3 /2 = 0.0053. (135) 
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FIG. 13: 

a) Diagrams, corresponding to the coefficient D\ in front of the term M 8 + M 8 in the Ginzburg-Landau expansion. The 
combinatoriai coefficient of the diagrams is 1/8. All the loops are symmetric under the replacement Q — > — Q, which is 
equivalent to the replacement u — ► v. Therefore all such diagrams are equal. 

b) Two different diagrams, which form the coefficient D2 in front of the term |w| 2 |it| 6 . The combinatorial coefficient of both is 
equal to 1. 

The diagrams c), d), e) form the coefficient D3 in front of the term |it| 4 |u| 4 : 

c) two loops with equal values but different combinatorial coefficients 1 and 1/2 correspondingly. Under the replacement u — * v 
both diagrams do not change. The complex conjugate of them simply reverse the encircling along the loop and such diagrams 
are taken into account in the combinatorial coefficients. 

d) the combinatorial coefficient of these both two diagrams is equal to 1, if we keep in mind the complex conjugate of them. 
The two loops transfer one into another under the replacement u — ► v. 

e) a diagram with a combinatorial coefficient equal to 1. 
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e) 



FIG. 14: Diagrams a), b) are of the sixth order in u: 

a) the Cooper loop for the third harmonic, which corresponds to the coefficient h in the GL expansion; it has a combinatorial 
coefficient equal to 1/2. 

b) diagrams, corresponding to the term 12(u-}qUq 2 v~q + v-zqv*_q 2 uq + h.c). The combinatorial coefficient of all four is equal 
to 1/2. 

Diagrams c), d), e) are of the eighth order in u: 

c) diagrams with a combinatorial coefficient of all four of them equal to 1. The upper two correspond to the term h{\uiQ | 2 \uq | 2 + 
\ V -3Q | 2 \ V -Q | 2 )i an d the lower two correspond to the term U(\u3Q \ 2 \v~q \ 2 + «-3Q | 2 \uq | 2 ). 

d) diagrams, corresponding to the term h(u3QV-zQUqV*_Q + h.c). The lower loop is a complex conjugate of the upper loop. 
The combinatorial coefficient is equal to 1. 

e) diagrams, corresponding to the terms IqU3qUq 2 V-q \uq | 2 and (I7 + Is)u3qUq 2 V-q\v-q\ 2 . There are even more loops which 
are complex conjugate to the ones shown on the picture, as well as loops obtained under the replacement u — > v. The values 
of such loops are the same, therefore we do not show them. 
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